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MODEL BASED FAULT DETECTION AND 
DIAGNOSIS METHODOLOGY FOR HVAC SUBSYSTEMS 

Abstract 

A fault detection system for an HVAC system including sensors for 
measuring the performance of a condenser, a compressor, an evaporator, and a chiller. 
The fault detection system is provided with a thermodynamic preprocessor for 
calculating characteristic quantities (CQ's) from a plurality of measured inputs. Also 
provided is a base-case lookup table for storing plural first sets of CQ values 
generated by the thermodynamic preprocessor in an initial period in which fault- free 
operation of the HVAC system is assumed, a first set of CQ values being generated 
for each of plural different measured input values. An interpolator is provided for 
interpolating a set of base-case CQ values from the first sets of CQ values stored in 
the base-case lookup table for a given set of measured inputs. A fault-detector detects 
a fault when a difference between actual CQ values and base-case CQ values exceeds 
a predetermined threshold value for at least one of the CQ values. A fault classifier 
classifies a detected fault based on which ones of the actual CQ values exceed the 
interpolated CQ values. 
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MODEL BASED FAULT DETECTION AND 
DIAGNOSIS METHODOLOGY FOR HVAC SUBSYSTEMS 



FIELD OF THE INVENTION 

1 The present invention relates to a method for detecting and diagnosing faults 

2 in HVAC subsystems. More specifically, the method of the present invention detects and 

3 diagnoses faults using a combination of a physical (thermodynamic) model and a neural 

4 network. 

5 BACKGROUND OF THE INVENTION 

6 Vapor-compression systems account for a very significant portion of energy 

7 consumption in the industrial and commercial sectors. In large office buildings for example, 

8 it is estimated that 1 0% to 25% of the total electric consumption can be attributed to cooling 

9 systems alone. Furthermore, these percentages can be much greater if a chiller subsystem 
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1 is operating at low performance levels due to the presence of faults . The extensive research 

2 in Fault Detection and Diagnosis (FDD) in HVAC systems has been motivated by concerns 

3 that range from: the need to reduce power consumption and consequently energy costs; to 

4 improve comfort levels in buildings; to reduce wear on various HVAC equipment; to reduce 

5 the magnitude of greenhouse emissions; and the need to support building operators in 

6 decision-making for building optimization. 

7 Conventional Energy Management and Control Systems (EMCS) are limited 

8 in that they do not provide the operator with the tools to diagnose faults in real-time. If 

9 effective FDD tools are designed and implemented in EMC systems, then the detection and 

1 0 diagnosis of various types of faults can be done automatically. 

1 1 A fault, in the context of HVAC applications, is defined as an unsatisfactory 

12 or unacceptable condition in the operation of a system or subsystem. A condition is 

1 3 unacceptable if it is a failure or if it causes one directly or indirectly through a series of other 

14 faults. There are various types of faults, some of which are more difficult to detect and 

15 diagnose than others. The three basic fault categories, ranked in order of severity, are 

16 degradations, malfunctions and hard failures. 

17 Degradation faults occur at a gradual rate but progressively worsen over a 

1 8 period of time. Overall system performance is not drastically changed until the degradation 

19 has matured beyond a critical level. Examples of system performance degradation faults are 

20 the fouling of condenser or evaporator tubes, or the clogging of filters. Notably, dirt and 

2 1 grime builds up on the heat transfer surfaces of the heat exchangers over time, which in turn 

2 
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1 modifies the heat transfer coefficient of these devices. These faults gradually result in 

2 pressure drops, reduced flow rates, or higher temperature differences between the two fluids. 

3 Malfunction faults have more immediately noticeable end effects than 

4 degradation faults, and tend to justify more immediate service. An example of this category 

5 of faults is valve leakage, sensor errors, controller breakdowns and damper stoppages are 

6 other examples. These faults usually result in failure to maintain the desired set points in 

7 temperature, pressure or flow rates. 

8 Seized compressors and broken fan belts are examples of "hard" or "complete" 

9 failure faults. In such circumstances the impact of the fault is severe, and justifies the most 

10 immediate service. Building operators must often shut down entire subsystems in order to 

1 1 carry out the necessary repair. 

12 Fault Detection and Diagnosis (FDD) is the initial step towards taking 

13 corrective or preventative measures in an HVAC system. Detection and diagnosis sometimes 

14 overlap and other times are treated separately. Fault detection involves the determination 

15 that a fault truly exists based on an observable quantity exceeding some predetermined 

16 threshold or criterion. The choice of this criterion or threshold is a tradeoff between the 

1 7 sensitivity of the FDD scheme and the likelihood of sounding a false alarm. Fault diagnosis 

18 is the subsequent step, that isolates the cause of the fault. Sometimes this diagnostic step 

1 9 does not immediately locate the root cause of the fault initially detected but may involve a 

20 series of steps that eventually converges on the cause at some later point in the diagnosis. 
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1 Once the fault has been diagnosed, fault evaluation is performed to assess 

2 whether the impact of the fault sufficient to require immediate service. In general, service 

3 should be performed whenever: (1) Comfort cannot be maintained; (2) Equipment or 

4 personal safety is compromised; (3) Environmental damage is done (e.g. refrigerant leakage); 

5 and (4) Service expense is justified by reduced energy costs. 

6 There are several approaches to fault detection. Most, if not all, involve 

7 observation of the differences between an "actual" quantity of the system (or component) 

8 during its operation and some "predicted" value of this same quantity. The quantity under 

9 scrutiny is sometimes the "raw" measured variables (e.g., temperatures, pressures or humidity 

1 0 ratios that have not been pre-processed) or it may be the result of some preprocessing that 

11 transforms the measured variables collectively (in specific subsets) or individually. Some 

12 examples of the quantities used for comparison in fault detection arc: Measured equipment 

13 performance vs. Model-based prediction of equipment performance; Measured equipment 

14 performance vs. Common sense expectations of acceptable performance; Measured 

1 5 thermodynamic states (e.g., temperatures, pressures, humidity ratios, etc.) vs. Model-based 

16 predicted thermodynamic states; Calculated (estimated) physical parameters for "current" 

17 performance vs. Calculated (estimated) physical parameters for "baseline" performance. 

1 8 Similarly, there are different approaches to fault diagnosis. The result of the 

19 diagnosis is not binary as in fault detection (i.e., fault or no-fault) but involves a selection 

20 from a range of different possibilities. One such fault diagnosis approach is a rule-based 

21 diagnosis, which involves a set of rules such as directional change in temperature deviations 

■ 4 
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1 that may be unique to a given fault. Another approach involves comparison of physical 

2 parameters for "current" operation and "normal" operation. For example, the conductance- 

3 area product, UA of a condenser can be calculated from entering and leaving temperatures, 

4 and may be used to diagnose condenser fouling. Yet another approach uses a pattern 

5 recognition technique applied to "current" residuals and a matrix of "expected" residual 

6 changes associated with each fault. 

7 Different kinds of models have been used in FDD methods. However, 

8 conventional models are unsatisfactory as they have limited general applicability. Notably, 

9 conventional models require extensive, time consuming tuning in order to tailor the model 

10 to a specific implementation. As a result, such models are not readily portable to different 

1 1 HVAC implementations. 

12 Accordingly, one object of the present invention is to provide an improved 

1 3 model-based FDD methodology that uses a combination of a physical model and a black-box 

14 model which has general applicability to HVAC systems. 

15 A further object of the present invention in this research is to develop an 

1 6 improved model-based FDD methodology for determining when and where problems occur 

17 in a centrifugal chiller system, and that is suitable for online implementation using day-to- 

18 day operating field data. 

1 9 Another object is to provide an improved model which compares characteristic 

20 quantities for everyday operating conditions to a predetermined baseline of a chiller's normal 
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1 (fault-free) operation, and uses statistical hypotheses to determine if a significant deviation 

2 from this normal operation has occurred. 

3 Yet another object is to provide a FDD method which accounts for sensor 

4 errors to provide a more robust FDD methodology. 

5 SUMMARY OF THE INVENTION 

6 The above-listed objects are met or exceeded by the present fault detection and 

7 diagnosis system which includes sensors for measuring the performance of a condenser, a 

8 compressor, an evaporator, which are all components of the chiller subsystem. The fault 

9 detection system further includes a thermodynamic preprocessor, a base-case lookup table, 

10 an interpolator, a fault-detector and a fault classifier. 

11 The thermodynamic preprocessor, which is preferably a software routine 

1 2 executed on a computer or the like, calculates characteristic quantities (CQ's) from a plurality 

13 of measured inputs. The base-case lookup table stores plural first sets of CQ values 

14 generated by the thermodynamic preprocessor in an initial period in which fault-free 

1 5 operation of the HVAC system is assumed, a first set of CQ values being generated for each 

1 6 of plural different measured input values. 

1 7 The interpolator generates a set of interpolates base-case CQ values for a given 

1 8 set of measured inputs from among the first sets of CQ values stored in the base-case lookup 

19 table that is stored in memory of the computer. The fault-detector, which is preferably a 

20 software routine executed on a computer or the like, detects a fault when a difference 

6 
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1 between actual CQ values, calculated by the thermodynamic preprocessor for a given set of 

2 measured inputs, and the interpolated set of base-case CQ values for the given set of 

3 measured inputs exceeds a predetermined threshold value for at least one of the CQ values. 

4 The fault classifier, which is preferably a software routine executed on a computer or the 

5 like, classifies a detected fault based on which ones of the actual CQ values exceed the 

6 interpolated CQ values. 

7 BRIEF DESCRIPTION OF THE DRAWINGS 

8 FIGURE. 1 is block diagram showing the overall structure of the model-based 

9 fault detection and diagnosis system of the present invention; 

10 FIG. 2 is a block diagram showing the inputs and some of the outputs of the 

1 1 thermodynamic preprocessor of the present invention; 

12 FIG. 3 is an EES (Engineering Equation Solver) procedure used to determine 

1 3 the logarithmic mean temperature difference between the evaporator saturation temperature 

14 (T_evap), the chilled water supply temperature (T_chws), and the chilled water return 

15 temperature (T_chwr); 

16 FIG. 4 is an EES procedure used to determine the logarithmic mean 

17 temperature difference between the condenser saturation temperature (T_cond), the 

18 condenser water supply temperature (T_cws) and the condenser water return temperature 

19 (T_cwr) 
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1 FIG. 5 depicts equation sets for the condenser side input data used by the EES 

2 procedure of FIG. 4; 

3 FIG. 6 are EES program statements for inputting data used in the equations of 

4 FIGs. 3 and 8; 

5 FIG. 7 are EES program statements for inputting data used in the equations of 

6 FIG. 9 and FIG. 10; 

7 FIG. 8 are equations used to determine evaporator related characteristic 

8 quantities; 

9 FIG. 9 are equations used to determine compressor related characteristic 
10 quantities; 

H FIG. 10 are equations used to determine condenser related characteristic 

12 quantities; 

13 FIG. 1 1 is an expansion valve equation used by the EES code in FIG.8; 

14 FIG. 12 is a list of thermodynamic properties and conversion constants used 

15 by the EES code in FIGs. 8-10; 

16 FIG. 13 is a schematic of a chiller showing the location of each of the sensors 

17 used to obtain the measured inputs; 

1 8 FIG. 14 illustrates the structure of a single neuron unit; 

1 9 FIG. 1 5 shows the general architecture of a interconnected neural network with 

20 several neurons, links and layers; 
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1 FIGs. 1 6A-C are graphs illustrating how different smoothing parameters are 

2 used in smoothing a trained GRNN; 

3 FIG. 17 shows the neural network architecture of the GRNN algorithm of the 

4 present invention; 

5 FIG. 1 8 shows the sensitivity of UA when used to detect flow rate reductions; 

6 FIG. 19 shows the sensitivity of the entering/leaving water temperature 

7 difference when used to detect flow rate reductions; 

8 FIG. 20 shows the sensitivity of UA when used to fouling of the heat 

9 exchanger tubes; 

I o FIG. 2 1 shows the sensitivity of APPR when used to detect fouling of the heat 

I I exchanger tubes; 

1 2 FIG. 22 shows the sensitivity of rj isaampb . when used to detect compressor faults; 

13 FIG. 23 shows the sensitivity of Tj motor when detecting motor or transmission 

14 faults. 

1 5 FIG. 24 shows the sensitivity of UA C , and APPR C when detecting fouling faults; 

16 FIG. 25 shows the sensitivity of compressor isentropic efficiency when 

17 detecting compressor faults; 

1 g FIG. 26 shows the sensitivity of motor efficiency when detecting motor faults; 

19 and 

20 FIG. 27 is a chart summarizing the relationships illustrated in FIGs. 1 8-26. 
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1 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

2 The inventor of the present invention discovered that it is difficult to detect and 

3 diagnose a fault in a chiller subsystem from raw data collected from sensors. Notably, it is 

4 difficult to set a meaningful fault threshold which is applicable to different chill subsystem 

5 implementations, i.e., is not customized for a particular implementation. To overcome this 

6 problem, the inventor of the present invention utilize a series of characteristic quantities 

7 which are calculated from the raw data using well known thermodynamic principles. The 

8 use of such characteristic quantities eliminates the need to individually tune the system of the 

9 present invention for any given HVAC installation. 

10 Importantly, the inventor of the present invention discovered that some 

11 characteristic quantities (CQs) are more sensitive than others in detecting certain faults. 

12 Table I below, summarizes the outputs (characteristic quantities) and the inputs used to 

13 determine them. 

14 The sensitivity of each CQ value to a given fault was determined by simulating 

1 5 various faults in a controlled experiment. Each fault condition in addition to a base case (no- 

1 6 fault) condition was simulated. Table II below pairs particular fault conditions with potential 

1 7 CQs that may be used for their detection. The inventor of the present invention developed 

1 8 a system for detecting and diagnosing faults utilizing the relationship summarized in Table 

19 II. 
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TABLE I 







(Raw Data) || 


Out 


puts (Characteristic Quantities) | 


1 


— 


chilled water supply 
temp. (°F) 


APPR e 


evaporator heat exchanger approach 
CD 


2 





chilled water return 
temp. (°F) 


APPR C 


condenser heat exchanger approach 
(°F) 


3 


T cws 


condenser water 
supply temp. (°F) 


CHWTD 


chilled water temp, difference (°F) 


4 


m chw 


chilled water flow rate 
(gpm) 


CWTD 


condenser water temp, difference (°F) 


5 


m cw 


condenser water flow 
rate (gpm) 


UA C 


condenser conductance-area product 
(Btu/hr-°F) 


6 


T evap 


evaporator saturation 
temp. (°F) 


UAe 


evaporator conductance-area product 
(Btu/hr-°F) 


7 


Tcond 


condenser saturation 
temp. (°F) 


"Hiseatropic 


isentropic heat transfer coefficient 
(Btu/hr-ft 2 -^) 


8 


T 2 


compressor discharge 
temp. (°F) 


*1 motor 


motor/transmission heat transfer 
coefficient (Btu/hr-ft 2 -°F) 


9 


P 


compressor power 


COP ovcrall 


overall coefficient of performance 


10 






cop £U111P 


compressor coefficient of performance 



TABLE II 



13 


Fault Type 


Sensitive CQs 


14 


(1) Condenser water flow reduction 


UA C , APPR C , CWTD 


15 


(2) Chilled water flow reduction 


UA C , APPR C , CHWTD 


16 


(3) Evaporator tube fouling 


UA e and APPR e 


17 


(4) Condenser tube fouling 


UA C and APPR C 


18 


(5) Compressor internal faults 


P, COP comD and % en , rowc 


19 


(6) Motor/Transmission faults 


P,COP overa „and7 MO , or 



FIG. 1 is a block diagram of the overall structure of the model-based fault 
detection and diagnosis methodology of the present invention. The fundamental components 
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1 along with their links to other components and flow of information for fault detection and 

2 diagnosis are depicted. 

3 There are two sources of data used in this methodology, operating CQ's 

4 calculated on-the-fly, and base-case CQ's which are interpolated from data stored in a look- 

5 up table. As shown in FIG. 1, measured inputs (and forcing functions), such as the chilled 

6 water supply temperature of the chiller and cooler water are fed to both the thermodynamic 

7 preprocessor 1 00 and the interpolator 1 02. The thermodynamic preprocessor 1 00 uses the 

8 measured inputs to generate the characteristic quantities (CQs). The operating CQ's are 

9 compared with base-case CQ's which are interpolated from empirical data stored in a base 

10 case lookup table 104. 

11 The look-up table 104 stores values of CQ's that were generated during a 

12 period in which fault-free operation of the chiller is assumed. The interpolator 102 

13 interpolates normal (base-case) CQ's from the data stored in the look-up table 104. The 

14 process for obtaining the operating and base case CQ values will be explained in detail 

15 below. 

16 The streams of operating and base-case CQ data are compared and errors 

1 7 (residuals) are computed. The residuals are fed to the Fault Classifier section 1 08 where the 

1 8 appropriate statistical analyses are employed to determine the significance of these errors. 

19 If they are indeed statistically significant, a fault detection alarm is produced. Otherwise, 

20 monitoring continues with the assumption that the chiller process is operating normally and 

21 that no faults are present. 

12 
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1 If faults are detected, the diagnostic classifier determines the type of fault and 

2 its likely location within the chiller. The CQs used in the methodology of the present 

3 invention were specifically selected to make isolating the fault location easier. As noted 

4 above, the inventor of the present invention discovered that certain CQs are sensitive to some 

5 faults and not to others. Thus, the careful selection of CQs makes it possible to distinguish 

6 between fault types. 

7 Once faults are detected and diagnosed, they are subsequently evaluated to 

8 determine the most suitable action to take. Fault Evaluation deals with comfort, safety, 

9 environment and economic issues. Action could be taken whenever: (1) comfort cannot be 

10 maintained, (2) equipment or personal safety is compromised, (3) environmental damage 

1 1 occurs (e.g. refrigerant leakage), or (4) reduced energy costs justify the service expense. 

12 In order to conduct on-line FDD using day-to-day operating data, it is desirable 

13 to have a model that can generate CQ's without the need for any manual intervention. The 

1 4 thermodynamic preprocessor 1 00 requires only eight measurement variables to generate the 

15 necessary CQ's immediately. See FIG. 2. Sets of CQ's are generated at two different time 

16 periods. As will be explained below, plural first sets of CQ variables (base-case CQ's), 

17 representative of normal (fault-free) system operating conditions, are generated and stored 

18 in the look-up table 1 04 during a preliminary period when a sufficiently wide range of loads 

19 is covered. A second set of CQ data (operating CQ's) is generated "on-the-fiy", i.e., 

20 generated contemporaneously with the receipt of the measured inputs to reflect the existing 

13 
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1 state of the system, after the preliminary period has expired, and is compared with the first 

2 set of CQ variables to determine if faults have occurred. 

3 Both sets of CQ variables are generated in the same manner, the only 

4 difference being the time period in which each set is generated. Namely, the plural sets of 

5 base-case CQ's are generated during a preliminary period in which fault-free operation of the 

6 chiller is carefully controlled. The operating set of CQ's are generated on-the-fly after the 

7 base-case CQ sets have been stored in the look-up table 104. 

8 The thermodynamic preprocessor 100 generates the operating CQ's using a 

9 thermodynamic data reduction program which was written in EES, an Engineering Equation 

10 Solver Software developed by S.A. Klein and F.L. Alvarado, F-Chart Software, Inc. The 

1 1 EES program is particularly suited for solving thermodynamic equations; however, one of 

12 ordinary skill in the art will appreciate that other programming languages or software 

13 packages may also be used. 

14 The thermodynamic data reduction program consists of two procedures (FIGs. 

15 3 and 4) and equation sets for the condenser side input data (FIG. 5), evaporator side input 

16 data (FIG. 6), compressor input data (FIG. 7), evaporator equations (FIG. 8), compressor 

17 equations (FIG. 9), condenser equations (FIG. 10), expansion valve equations (FIG. ll)and 

18 forthethermodynamicproperties and conversion constants (FIG. 12). It should be noted that 

19 the thermodynamic data reduction program uses general thermodynamic principles, and is 

20 not specifically timed for a specific chiller subsystem. For this reason, the CQ's have general 

2 1 applicability to all chiller subsystems . 

14 
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1 A procedure statement heads each procedure (FIGs. 3 and 4), and identifies the 

2 name and arguments of the procedure. Thus, the statement "Procedure LMTD_e(T_chws, 

3 T_chwr, T evap: LMTDe, T_we, marge)" defines a procedure LMTD_e and a series of 

4 arguments (inputs and outputs). 

5 Table I, above, lists the input/outputs of the thermodynamic data reduction 

6 program. Moreover, FIG. 13 is a schematic of a chiller showing the location of each of the 

7 sensors used to obtain the measured inputs. 

8 The procedures (FIGs. 3 and 4) utilize function calls in which the function 

9 max(x,y) returns the larger of the two arguments (x, y). Thus, max(0, -2) returns "0". 

10 Correspondingly, the function ln(x) returns the natural log of argument x. 

11 The procedure LMTD_e (FIG. 3) is used to determine the logarithmic mean 

1 2 temperature difference between the evaporator saturation temperature (T_evap), the chilled 

1 3 water supply temperature (T_chws), and the chiller water return temperature (T_chwr). The 

14 built-in EES error procedure halts the calculations if T_evap is greater than T_chws or if 

15 T_chws is greater than T_chwr. 

16 The procedure, LMTD c (FIG. 4) is used to determine the logarithmic mean 

17 temperature difference between the condenser saturation temperature (T_cond), the 

1 8 condenser water supply temperature (T_cws) and the condenser water return temperature 

19 (T_cwr). The built-in EES error procedure halts the calculations if T_cond is less than 

20 T_cwr or if T_cwr is less than T_cws. 
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1 In developing the base case look-up table, measurement data for a wide variety 

2 of base-case operating conditions is compiled in a table which is accessed by the 

3 thermodynamic data reduction program. As will be appreciated by one of ordinary skill in 

4 the art, it is advantageous to collect measurement data for as wide a variety of operating 

5 conditions as possible in order to improve the accuracy of the interpolated CQ values. For 

6 example, it is desirable to calculate base-case CQ values for various times of day for both 

7 workdays and non-workdays, different weather conditions (ambient temperature, sunlight, 

8 winds, humidity), etc to fully reflect the wide range of heating/cooling loads encountered by 

9 the system. 

I o In FIGs. 5-7, the statement lookup( x, y) "F" is used to access entry "x" in the 

I I lookup table and assign the returned value to variable "y". The units for each variable 

12 accessed from the lookup table is provided in the comments following each lookup 

13 statement. Thus, in the above example, the units are in degrees Fahrenheit. 

1 4 As noted above, a first set of CQ variables, representative of normal (fault-free) 

15 system operating conditions, is generated during a preliminary period when a sufficiently 

1 6 wide range of loads is covered. It is assumed that sufficient data are collected and that there 

17 are no major unforeseen systematic changes during data collection. 

18 Fault detection is performed by comparing each respective operating CQ 

1 9 variable with a corresponding base-case CQ variable. Each set of CQ variables stored in the 

20 lookup table 104 includes at least eight (of the nine) input and ten (of the output variables 

21 listed "in Table I. 

16 
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1 As will be appreciated by one of ordinary skill in the art, there are many 

2 different combinations of input variables, whereas the look-up table only contains a relatively 

3 small subset of the possible combinations. Consequently, it is unlikely that a given set of 

4 eight measured inputs acquired on-the-fly will directly match one of the sets of input 

5 variables stored in the look-up table 1 04. Accordingly, the interpolator 1 02 determines a set 

6 of base-case CQ variables through a process of interpolation. 

7 The inventor of the present invention investigated a number of different 

8 interpolation schemes. Linear interpolation and regression against power form relations were 

9 found to be inaccurate and not general. According to a preferred embodiment, Neural 

10 Network architecture was found to be well suited to the interpolation. However, one of 

1 1 ordinary skill in the art will appreciate that other interpolation methods may be substituted. 

12 Neuralnetworkmodelsareinspiredbythehumanthoughtprocesses. Over 100 

13 billion biological neurons are present in the human brain. The connections between these 

14 neurons are called synapses and when the brain learns their strength is modified. 

15 Analogously, artificial neural networks contain artificial neurons that are 

1 6 connected via one-way information conduits called links. Weights are associated with these 

17 links that control the magnitude of the input signal entering the artificial neuron. These link 

1 8 weights simulate the physical and neuro-chemical characteristics of the biological synapse. 

1 9 Each artificial neuron in the network functions by first summing its scaled inputs and then 

20 applying a non-linear function to this sum to generate an output. FIG. 14 illustrates the 

17 



CA. 02344908 2001-04-23 

1 simple structure of a single neuron unit and FIG. 15 shows the general architecture of a 

2 interconnected network with several neurons, links and layers. 

3 A layer is a grouping of neurons and are of three basic types. The input layer 

4 receives data from outside of the network, the output layer holds the network's final 

5 computational results and any layers in between these two are called hidden layers. 

6 A basic feature of neural networks is that they learn by example and are never 

7 programmed. Therefore, in order to train a network, it must be presented with a data sample 

8 thatcontainstheinputsandtheircorrespondingoutputs. By an iterative process, the network 

9 gradually learns the input-output relationship andthen, depending on the application, is used 

10 for prediction, correction or pattern recognition. Once trained, they can resolve numeric 

1 1 problems that are otherwise solved by conventional regression analysis. In the network, the 

1 2 inputs are equivalent to the independent variables used in regression, the dependent variables 

13 become the outputs and the observations are the sample patterns used for training the 

14 network. 

! 5 Neural networks fall into the two basic categories of: supervised networks and 

16 unsupervised networks. Supervised networks are those that make predictions, classify 

17 patterns or make decisions based on a set of learned inputs and outputs. Some examples of 

18 t hesearetheProbabilisticNeu ra lNetworks(PNN),theGeneralRegressionNeuralNetworks 

19 (GRNN) and Back Propagation Networks. These kinds of networks are trained to make 

20 predictions, classifications, or decisions when given a large enough number of accurate 

21 classifications or predictions from which it can learn. 

18 
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Unsupervised networks, without being taught how to categorize, are able to 
classify a given number of input and output patterns into a specified number of categories. 
They do this by clustering the training patterns using their proximity in « dimensional space 
where n is the number of inputs. The network usually clusters the data into the maximum 
number of categories it is presented with. Kohonen Networks are examples of this network 
type. 

The GRNN is one type of neural network that is well suited to interpolation. 
It is based on the estimation of a probability density function of a vector random variable, 
X, and a scalar random variable, Y. If the joint probability density function of these variables 
are both known then the conditional probability density function and the expected value can 
be computed. The estimated value of Ffor a given X is presented in the following general 
regression equation: 

\Yf{X,Y)dy 

E[Y\X] = ^~ (1) 

\f(X,Y)dy 

where 

l{Y\X] = conditional mean of Y on X 

f(X, Y) = known joint continuous probability density function 
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The probability density function,/(X, Y), is estimated from sample observations 
of X and Y when it is unknown. For a non-parametric estimate of f(X,Y), the Parzen 
estimation f(X, Y), is used by the GRNN. It is defined by the following equation for the 
observed sample observations, X- t and Y t of the vector X and scalar Y: 



where 



f'(X,Y)=- 



1 ~ztf*f> 



(2) 



/ A ,=exp 



(x-xfjx-x,) 

2o> 



(3) 



/ r =exp 



2o> 



(4) 



and 

n = the number of sample observations 
p = the dimension of the X vector 
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a - the standard deviation (or smoothing parameter) 

An estimate for the desired mean of Fat any given Xis derived in Equation (5) 
by combining Equations (1) and (2) and performing the integration after first interchanging 
the integration and summation operation. 




(5) 



where, the scalar function is given by: 

D?={X-X$(X-X) (6) 

The main algorithm of the GRNN model is expressed by Equations (5) and (6). 

The estimate Y(X) is a weighted average of all the observed samples, Y> where each sample 

is weighted in an exponential manner according to the Euclidean distance, D„ from each^.. 
This appropriate weighting is explained by the inversely proportional relationship between 

theexpression ex pf-£j] and A- That is, as £>,. increases, e ^-£j) decreases and vice- 



u" An optimum value for the smoothing parameter, a, is determined using iterative 
21 
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1 genetic adaptive methods. Iterative methods such as the Holdout method and the Wiggle 

2 method are well known by those of ordinary skill in the art, and do not form part of the 

3 present invention. Larger values of a improve smoothness of the regression surface. 

4 Notably, a must be greater than 0 and can usually range from .01 to 1 with good results. 

5 FIGs. 16A-C illustrates how different smoothing parameters are used in 

6 smoothing a trained GRNN. In this example, one input value is used to predict one output 

7 value. As shown in the graph, of the forty input and output patterns only input patterns 1 0, 

8 20, and 30 produce an output value of 1 . 

9 FIG. 17 represents the neural network architecture of the GRNN algorithm of 

10 Equations (5) and (6). The Euclidean distance, £>,., is computed by the links between the 

11 input layer and the first hidden layer. Based on observed samples, X t , and smoothing 

12 parameter, a, the expression, 



13 




14 is computed. A node in the second hidden layer takes the sum of the exponential values of 

1 5 all samples. In other nodes of this same layer, the products of the exponential values and the 

1 6 corresponding observed Y- t for each sample observation are computed. The node in the third 

17 hidden layer computes the sum of all these product values, which is then supplied to the 

1 8 output node where the ratio between it and the previous sum is calculated. 
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1 The suitability of the GRNN, however, is attributed to several important 

2 features and makes it convenient for online implementation. In the GRNN network only a 

3 single parameter is estimated. Unlike other networks, a once through, non-iterative training 

4 process with a highly parallel structure is involved. Compared to conventional regression 

5 analysis, the specification of the underlying regression function, bounds of the independent 

6 variables, initial convergence values and convergence criteria are not required beforehand. 

7 Additionally, the algorithm provides smooth transitions from one observed value to another 

8 even with sparse and noisy data in a multidimensional measurement space and can be used 

9 for any regression problem where an assumption of linearity is not justified. 

10 Since the thermodynamic and the GRNN models involve functions and 

1 1 variables in a multi-dimensional space, it is important to know what the truly independent 

12 measurement variables are. Here, the term independent refers to those variables that arc 

13 forcing functions of the chiller process. Information about variable independence is 

14 necessary for understanding the available data, achieving greater insight to data 

15 manipulation, accurately assigning inputs and outputs, and successfully training the GRNN 

16 model to be able to predict normal system behavior. For the chiller process, only five 

17 measurement variables may be considered independent This number is determined from 
1 s studying the various equations that define the chiller process more closely in the discussion 

19 that follows. 

20 An overall energy balance of the entire refrigeration system is: 
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j Power =Q atp -Q ami (?) 

3 4^=»V(^-^) = "V^,a«f (9) 

4 The enthalpy difference, Ah refcond , between the entering and leaving refrigerant 

5 in the condenser is a function of the condenser saturation temperature, T cond : 

6 Ah nffi0lxl =f(T C0Kl ) ^ 

7 Q*=<£JT*,-T+.) (11) 

8 a^=^(^4-^)=<^^ ■ (12) 

9 Here, the enthalpy difference, At] rehevop , between the entering and leaving 
1 0 refrigerant in the evaporator is a function of the evaporator saturation temperature, T^: 

1 2 In addition to the above relations, the following relations for the chiller Power 

13 and Capacity, are available from chiller manufacturers' experiments: 
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1 Pomt^KT^TcoJ (14) 

2 Qe^^fi^^cond) (15) 

3 Equations (7) to (15) folly define the thermodynamic data reduction model. 

4 An inventory of the number of equations and the number of unknowns yield nine and 

5 fourteen, respectively. Thus the difference between these two numbers tells us that five 

6 variables must be independently measured for these model equations to work. 

7 The inventor of the present invention selected the following five independent 

8 measurement points condenser supply temperature, T cws , the chilled water supply 

9 temperature, F cAhij , the chilled water return temperature, T chwn the condenser water flow rate, 

10 and the chilled water flow rate, m^. However, the selection of other independent 

1 1 measurement points is contemplated and falls within the scope of the present invention. 

12 The condenser water supply temperature, T cws and the chilled water return 

1 3 temperature, T chwr are the entering water temperatures to the chiller subsystem. The chiller 

14 subsystem has no direct effect on the magnitudes of these two water streams until after they 

1 5 cross the system boundary. Therefore, the temperatures, T chwr and are independent. They 

1 6 represent only a part of the building load and heat rejection, respectively. The flow rates, 

17 fn^ and are additional required independent variables. The final measurement 

25 
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1 variable, the chilled water supply temperature, T chw , is termed independent because its 

2 magnitude is associated with a control set-point 

3 The essence of the Fault Detection scheme relies on the determination that a 

4 fault truly exists based on the deviation of a Characteristic Quantity (CQ) exceeding some 

5 predetermined threshold or criterion. This deviation is measured from an established base 

6 case representing normal operating conditions. 

7 As described above, the CQ values are computed from measurement data 

8 gathered from sensors. Due to the many perturbations that exist, there is an inevitable chain 

9 of causes and effects that occur and form the basis for uncertainties in the measurement data. 

10 Some of the factors that may cause perturbations are power supply variance, hysteresis, 

1 1 precision and bias errors and general instability. These factors contribute to the accuracy 

12 ratings of the measurement sensors and therefore, are factors that contribute to the 

1 3 uncertainties in CQ deviations. These uncertainties in CQs determine the threshold and thus 

14 the sensitivity for fault detection. A tighter uncertainty region corresponds to more sensitive 

15 detection. 

16 Numerous situations involve the computation of an important quantity, Y, 

17 instead of its direct measurement. If such a quantity is determined from N other directly 

1 8 measured quantities X { , X 2 ,..., X N , a functional relation between these quantities and Y may 

19 be expressed as Y= f{X x ,X 2 ,...,X N ). This function / not only expresses a physical law but 
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a measurement process. In fact, all the quantities that could contribute a significant 
uncertainty to the measurement results are accounted for by f. 

Furthermore, if ^represents the estimate of the output quantity, Y, and x„ x 2 , 
. . . , x N are input estimates for the values of the N input measurements, X„ X 2 ,..., X N , then Y 
can be expressed as Y=f(x {1 X li ... y X N '). In order to obtain the combined standard 
uncertainty, U C (Y) of the computation result Y (an estimated standard deviation of the 

result), a first-order Taylor series approximation of Y=f(X l ,X 2 ,...,X N ) is performed. This 
gives the following expression known as the Law of Propagation of Uncertainty: 



(16) 



= sensitivity coefficients 
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= standard uncertainty associated with the input estimate, x f 



= estimated covariance associated with*,- and Xj 



If the individual measurement inputs are assumed uncorrected and 
random, then the covariance term is zero and U C (J) is simplified to: 



condenser arc its approach, APPR C , and conductance-area product, UA C . Using the 
simplified form of the uncertainty propagation law in Equation (17), expressions can 
be derived for the uncertainties in UA C and APPR C . 

Equation (18) that defines the conductance-area product of the 
condenser, UA C , in terms of four independently measured variables and a fluid 
property, C^. This property is the specific heat of water and actually varies with 
temperature conditions but was assumed constant for this application since the range 
of variation is insignificant. 




(17) 



By manner of example, two CQs that are sensitive to tube fouling in the 




(18) 
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The UA C is a /(T^T^T^mJ. Equation (16) was used to express the 
most probable uncertainty in UA C as: 



u c (UA c > 



where, 



dJA r -fnC 



(19) 



(20) 



(21) 



(22) 



(23) 



and, KW^T^AT^) and ttjhj are the uncertainties of the measured 

condenser water supply temperature, condenser water return temperature, refrigerant 
condenser saturation temperature and the water mass flow rate, respectively. These 
uncertainties are equivalent to the accuracy of the measurement sensors. 
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The approach, APPR C is a /CC* T co J . The most probable uncertainty 
in APPR C may be expressed as: 



u(APPR c )= 



^4 m-^ 



(24) 



atppR c 



=-i 



(25) 



(26) 



For the threshold analysis, only precision errors were used. The 
following sensor precision errors, developed using time series analysis, were used to 
evaluate the effect of error on CQ: 

-Temperature sensor: ±0.3°F and ±0.6°F 

•Water flow rate sensor: ±5% 

•Power kW-meter: ±10% and ±20% 
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1 Fault detection thresholds are illustrated in FlGs. 18 - 24. Two types 

2 of horizontal lines are used in defining these thresholds. The no-fault line represents 

3 no flow rate reduction (or increase), no fouling and no inefficiencies in either the 

4 compressor or the motor. The critical line was determined by the following 3-step 

5 procedure: 

6 1 . Locate the points where the uncertainty boundaries intersect the no-fault 

7 line. 

8 2. Project these points vertically to the points that intersect the 0.0°F curve 

9 (shown as black dots; for ±0.3 °F, refer to FIG. 1 8). 

I o 3 Draw horizontal lines, parallel to the abscissa, through these intersection 

II points. 

12 These horizontal lines are the critical lines and they define a region 

1 3 (shown as hatched in FIG. 1 8) within which it is possible for faults to go undetected. 

14 Beyond these critical lines (or thresholds), however, faults may be detected with more 

1 5 confidence despite the errors present in the measurement sensors. 

1 6 For example, water flow rate abnormalities (reduction or increase) are 

17 detected most effectively by deviations in the conductance-area product and the 

1 8 entering/leaving water temperature difference of the given heat exchanger (evaporator 

19 or condenser). FIG. 18 shows the sensitivity of UA when used to detect flow rate 

20 reduction's. The heavy solid curve represents the UA deviation in the face of no errors 

21 in the temperature and flow sensors (i.e., 0.0°F/0% of full scale, gpm). The dashed 
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and light solid curves represent the UA deviation due to errors of ±0.3°F/±5% gpm 
and±0.6°F/±5% gpm, respectively. 

Here, the critical lines are drawn using the uncertainty boundaries 
defined by ±0.6°F/±5% gpm. Based on a ±0.6°F/±5% gpm precision of these sensors, 
it can be concluded that- UA deviations between -23.6% and +23.6% would not 
indicate a fault. Flow rate reductions greater than 41.0% and flow rate increases 
greater than 68.0% would be needed to indicate faults. 

FIG. 1 9 is similar to FIG. 1 8 in terms of the precision errors represented 
by the heavy, dashed and light curves. However, the characteristic quantity used here 
is the entering/leaving water temperature difference instead of the conductance-area 
product. It is apparent from the critical lines drawn in FIG. 19 that the thresholds for 
detecting flow rate reduction is lower than those depicted in FIG. 1 8 . In others words, 
this characteristic quantity has a greater fault detection sensitivity. 

Quantitatively, here, it is concluded that water temperature difference 
deviations beyond -14.3% and +14.3% are used to detect flow rate increases greater 
than 16.0% and flow rate reductions greater than 13.0%, respectively. 

Fouling describes a decrease in the effectiveness of the heat exchanger 
tubes, and is detected most effectively by deviations in the conductance-area product 
and the approach of the given device. FIG. 20 shows the sensitivity of UA when used 
to detect fouling. A fouling factor of 1 represents no fault and in this study, is the 
normalized equivalent of 0.00075 hr-fr7BTU. A factor greater than 1 is an increase 
32 
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1 in fouling and lesser than 1 is a decrease. The critical lines indicate, that UA 

2 deviations beyond -19.6% and +1 9.6% are used to detect fouling factors greater than 

3 1 .6 and lesser than 0.3, respectively. 

4 FIG. 21 shows the sensitivity of the Approach when used to detect 

5 fouling. The critical lines indicates thatAPPR deviations beyond -1 7.6% and +17.6% 

6 are used to detect fouling factors lesser than 0.7 and greater than 1.3, respectively. 

7 These numbers, compared to those for UA, suggest that the APPR has lower fouling 

8 detection thresholds. 

9 Some internal faults of the compressor are detected most effectively by 

10 deviations in its isentropic efficiency, h isenlrvpic . FIG. 22 shows the sensitivity of 

1 1 Vhentropic wnen used to detect faults - A compressor fault factor of one represents no 

12 fault. In contrast to the fouling factor, a compressor fault factor greater than one is 

13 an improvement of the fault condition and a factor lesser than one is a worsening of 

14 the situation. The critical lines indicate that 7j isenlroi>ic deviations beyond -2.1% and 

15 +2.1 % are used to detect fault factors lesser than 0.8 and greater than 1 .2 respectively. 

1 6 Deviations in the motor efficiency, 7j molor , are used to detect faults in the 

1 7 motor or transmission housing. FIG. 23 shows the sensitivity of 7? motor when detecting 

18 motor or transmission faults. The heavy solid curve represents the Tj motor deviation 

19 when the kW-meter is error free (i.e., 0% of full scale, kW). The dashed and light 

20 solid curves represent the Tj moior deviation due to errors of ± 1 0% kW and ±20% kW, 
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respectively. A motor fault factor of one represents no fault and in this study, is 
equivalentto 0.85. However, factors greater than one represent worsening conditions. 

Here, the critical lines are drawn using the uncertainty boundaries 
defined by the ±20% kW error. Based on this kW-meter precision, it is concluded 
that Tj motor deviations less than -20% is used to detect fault factors greater than 1.3. 

Equations (27)-(29) below are similar in form to the previous three 
relations for the condenser, and are used to model the chilled water reduction. Here, 
the chilled water flow rate reduction causes an increase in the return temperature and 
a decrease in the evaporator temperature for a given load and chilled water set point 
temperature. 

>Ww, = (27) 
Qe^ = C lM (T dnu . - T dn J (28) 

where, m^jvu , factor ae , NEPASS, NETUBE and D e are the chilled water flow due to 

the flow reduction fault, the percentage reduction in chilled water flow rate, the 
number of evaporator tube passes, the number of evaporator tubes and the diameter 
of the evaporator tubes, respectively. 
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1 The CQs UA e ,APPR e and CHWTD show the most significant deviations 

2 from the no-fault base-case. The CQ relationships here are similar to those used for 

3 the condenser water reduction fault. CHWTD is the most sensitive and its sensitivity 

4 is unaffected by load changes. 

5 Equations (30) and (31) are used to model condenser tube fouling. 

6 There are some conflicting interpretations of the fouling factor in the literature. Here, 

7 fouling is defined as the resistance to heat transfer due to the tube material as well as 

8 to any thin film build-up on its surfaces. Equation (30) is used to simulate an increase 

9 (or decrease) of the level of fouling in the condenser tubes. It is apparent from 
10 Equation (3 1) that as fouling increases, UA C decreases. 

H Kj^K-f^fc oo) 

12 UA C =- X ^ (31) 

+ +R 

1 3 where,factor fc is the fault modification factor used to simulate degrees of fouling in 

14 condenser (unity e normal) andR c/ault is the resistance to heat transfer associated with 

15 the condenser tube material, including the fouling factor for faulty conditions. 

1 g The UA C , and APPR C are CQs that show the more significant deviations 

17 from the rio-fault base-case. FIG. 24 shows how these CQs deviate from the base- 

1 8 case condition as the fault condition changes (i.e., as fouling increases or decreases). 
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1 Similar to the flow reduction faults, the different CQs behave differently. For 

2 example, at a base-case load of 1 500 tons, a 1 00% fouling increase by the condenser 

3 tubes corresponds to a 27% decrease and a 54% increase in UA C and APPR C , 

4 respectively. However, in contrast to condenser flow reduction faults, CWTD shows 

5 no significant change. This behavior is important in distinguishing the fouling fault 

6 from the flow reduction fault. 

7 At 500 tons, the same fouling condition corresponds to a 29% decrease 

8 and 63% increase in UA C and APPR C , respectively. It is apparent from these results 

9 that the more sensitive CQ is the APPR C and that varying loading conditions affect it 

10 more than they did UA C . 

1 1 Similar trends exist for a 50% fouling decrease condition. CQs that 

1 2 increased before decreased and vise- verse. Fouling decrease information is as useful 

13 as fouling increase information. Here, FDD may be used to check the integrity of 

1 4 scheduled maintenance, repair or cleaning of fouled tubes. 

1 5 Equations (32) and (33) are used to model evaporator tube fouling. 

16 (32) 



17 
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where, factor^, is the fault modification factor used to simulate degrees of fouling in 
the condenser (unity b normal) and R e and R eJau „ are the resistance to heat transfer 
associated with the evaporator tube material, including the fouling factor for normal 
and for faulty conditions, respectively. 

The UA e mdAPPR e are CQs that show the more significant deviations 
from the no-fault base-case. The CQ relationships here are similar to those used for 
the condenser tube fouling fault. APPR e is more sensitive and its sensitivity is slightly 
affected by load changes. 

The Weisner correlation, based on the measured performance of 
centrifugal compressors, is used for modeling compressor internal faults. The 
modification of the reference polytropic efficiency, rj re{er , in Equation (34) is used to 
simulate possible faulty conditions within the compressor. 

^tfU = Prefer factOf omp ( 34 ) 

where, factor comp and rj referfau!( , are the compressor fault modification factor, and the 
modified 7 re/er due to faulty conditions, respectively. 

The compressor isentropic efficiency, rj isenlropic is the CQ that shows the 
most significant deviations from the no-fault base-case. FIG. 25 illustrates a greater 
decrease in the isentropic efficiency as the compressor fault factor decreases. Based 
on the manner in which this fault is simulated, a decrease in the fault factor means a 
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1 worsening compressor internal fault condition. Varying the loading conditions does 

2 not appear to have a great impact on the deviation of the isentropic efficiency. 

3 Whenever there is a fault in the motor or transmission system of the 

4 compressor, the overall power draw (P) is expected to increase compared to the actual 

5 power used for the compression process. The modification factor in Equation (35) 

6 simulates this increase. Equation (35) defines the combined efficiency of the motor 

7 and transmission. An energy balance on the compressor establishes the relation 

8 between the states of the refrigerant at the suction and discharge conditions. 

9 P m =P '-factor^ (35) 

10 U= p • (36) 

Oh* 

11 where, P faulp factor motor and bTj^p are the compressor power consumption due to 

12 faulty conditions, the motor fault modification factor and the change in refrigerant 

13 enthalpy between compressor's suction and discharge states, respectively. 

14 The motor efficiency is the CQ that shows the most significant 

1 5 deviation from the no-fault base-case. FIG. 26 shows that there is a greater decrease 

16 in the motor efficiency as the motor fault factor increases (i.e., as the power draw 

17 increases for a given loading condition). From this simulation, the motor efficiency 

1 8 deviation (not the motor efficiency) is independent of loading conditions. That is, at 
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1 any loading condition, the percentage deviations due to an increased fault condition 

2 will show the same magnitudes. 

3 Refrigerant loss results from leakage in chiller subsystems, which has 

4 the net effect of reducing overall system performance. It is often confirmed by 

5 bubbles present in the refrigerant (as seen from the liquid line sight glass), frost 

6 formation, abnormal superheat and subcooling temperatures, and abnormal unit 

7 operating pressures. A temperature difference between the low pressure liquid line 

8 temperature (entering the evaporator) and the suction line temperature (entering the 

9 compressor) was computed and used as a characteristic quantity to flag the occurrence 

1 0 of refrigerant leaks. 

11 FIG. 27 summarizes the diagnostic patterns exhibited by the CQs 

12 discussed earlier. The symbols "0", "+" and "-" signify no significant change, an 

13 increase and a decrease compared to base-case values, respectively. In the Type 

14 column, "N", "F" and "N/A" stands for Normal, Fault and Not Applicable, 

15 respectively. Two normal operating conditions called "Load change" and "Heat- 

16 rejection change" are included. Respectively, these are quantities representing a 

17 change in the chiller load and a change in the amount of heat rejected by the 

1 8 condenser. They are modeled by varying the chilled water return and condenser water 

19 supply temperatures, respectively. These conditions are included to show that it is 

20 possible to distinguish them from fault conditions using the appropriate set of CQs. 
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1 Although, these CQs effectively isolate faults on the component level, 

2 the diagnostic patterns help to further distinguish between the various fault conditions 

3 that may exist in a particular component. For example, the increase and decrease in 

4 the condenser approach, APPR C and the condenser conductance-area product, UA C , 

5 respectively, would signify that there is a problem in the condenser and not in the 

6 other components of the chiller sub-system. However, this same CQ change partem 

7 is both exhibited by an increase in condenser tube fouling and a decrease in the 

8 condenser water flow rate. To distinguish between these two conditions, FIG. 27 

9 shows that the condenser water flow rate decrease is accompanied by an increase in 

10 the condenser water temperature difference, CWTD, whereas the condenser tube 

1 1 fouling is not. 

1 2 For most performance monitoring studies, the COP overall is a CQ that is 

13 often used to determine significant performance changes. The COP overa „ is included 

14 in FIG. 27 to illustrate that although it gives some indication of changes within a 

15 chiller sub-system, by itself, it is quite vague in isolating faulty components or in 

16 diagnosing faults within a given component. Therefore, it must be used in 

17 combination with the other CQs, mentioned here, for useful diagnoses to be made. 

18 According to one modification of the above-described FDD 

1 9 methodology, the use of a base-case lookup table is eliminated. Notably, the use of 

20 parallel fault detection and diagnosis in which CQ residuals are derived from 

21 comparing a particular chiller to other identical (similar model, capacity and run-time) 
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1 chillers within the same HVAC system eliminates the need for a lookup .table. In 

2 parallel FDD, the data from a designated chiller is considered to be the base-case, and 

3 faults are detected by comparing the data from another chiller with the data from the 

4 designated chiller. In so doing, faults in a particular chiller or chillers are determined 

5 over a predetermined time period. 

6 The parallel FDD methodology is substantially similar to the above 

7 described (serial) FDD methodology in all relevant respects. As will be appreciated 

8 by one of ordinary skill in the art, in parallel FDD the CQ's from the designated 

9 chiller and the other chiller(s) are determined on-the-fly, and the difference (residuals) 

10 between the two streams of CQ's is used to detect faults in the previously described 

11 manner. 

12 One of ordinary skill in the art will appreciate that both of the above- 

13 described embodiments of the present invention utilize CQ's and thus do not require 

14 modification or tuning of the thermodynamic model. Accordingly, the FDD system 

15 of the present invention may readily be included in an existing facility management 

16 system, provided that the appropriate sensors are already installed. Use of the first 

17 embodiment merely requires the collection and storing of base-case CQ's whereas 

18 base-case CQ's are calculated on the fly according to parallel FDD. 

19 Another benefit of the present invention relates to threshold used to 

20 detect the occurrence of faults. Both of the above-described embodiments detect of 

21 faults by comparing operating CQ's to base-case CQ's. A fault is detected when one 
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1 or more CQ's exceed a predetermined threshold which is determined in relation to the 

2 sensor accuracy. Accordingly, the false detection of faults is minimized. 

3 Yet another benefit of the claimed invention relates to the classification 

4 of faults. Notably, the inventor of the present invention discovered that certain CQ's 

5 are more sensitive than others to the occurrence of faults, and that faults may 

6 accurately be diagnosed by examination of which CQ's exceed the predetermined 

7 threshold. 

8 While various embodiments of the present invention have been shown 

9 and described, it should be understood that other modifications, substitutions and 

10 alternatives are apparent to one of ordinary skill in the art. Such modifications, 

1 1 substitutions and alternatives can be made without departing from the spirit and scope 

12 of the invention, which should be determined from the appended claims. 

13 Various features of the invention are set forth in the appended claims. 
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WHAT IS CLAIMED IS: 

1 . A fault detection system for an HVAC system including sensors for 
measuring the performance of one or more of a condenser, a compressor, an 
evaporator, and a chiller, said fault detection system comprising: 

a processing means for generating characteristic quantities (CQ's) from a 
plurality of measured inputs; 

means for storing plural CQ values generated by said processing means in an 
initial period in which fault-free operation of the HVAC system is assumed, CQ 
values being generated for each of plural different measured input values; 

means for producing interpolated base-case CQ values from said plural stored 
CQ values for a given set of measured inputs; 

means for detecting a fault when a difference between actual CQ values, 
calculated by said processing means using said given set of measured inputs, and said 
interpolated base-case CQ values varies from a predetermined threshold value; 

means for classifying a detected fault based on which ones of said actual CQ 
values varies from said interpolated CQ values. 

2. The fault detection system according to claim 1 wherein said 
processing means uses measured inputs selected from the group comprising: chilled 
water supply temperature; chilled water return temperature; condenser water supply 
temperature; chilled water flowrate; condenser water flow rate; evaporator saturation 
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temperature; condenser saturation temperature; and compressor discharge 
temperature. 

3. The fault detection system according to claim 1 wherein said 
interpolator is a General Regression Neural Network. 

4. The fault detection system according to claim 1 wherein said 
interpolator is a Probabilistic Neural Network. 

5. The fault detection system according to claim 1 wherein said 
interpolator is a Back Propagation Network. 

6. The fault detection system according to claim 1 wherein said CQ's are 
selected from the group comprising: evaporator heat exchanger approach; condenser 
heat exchanger approach; chilled water temperature difference; condenser water 
temperature; condenser conductance-area product; condenser conductance-area 
product; isentropic efficiency; motor/transmission efficiency; overall coefficient of 
performance; and compressor coefficient of performance. 

7. ' The fault detection system according to claim 6 wherein condenser 
water flow rate abnormalities are diagnosed when one of said condenser conductance- 
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area product, said condenser heat exchanger approach, and said condenser water 
temperature difference exceeds said predetermined threshold value for said CQ. 

8. The fault detection system according to claim 6 wherein chilled water 
flow rate abnormalities are diagnosed when one of said evaporator conductance-area 
product, said evaporator heat exchanger approach, and said chilled water temperature 
difference exceeds said predetermined threshold value for said CQ. 

9. The fault detection system according to claim 6 wherein chilled water 
flow rate abnormalities are diagnosed when one of said evaporator conductance-area 
product, said evaporator heat exchanger approach, and said chilled water temperature 
difference exceeds said predetermined threshold value for said CQ. 

10. The fault detection system according to claim 6 wherein evaporator 
tube fouling is diagnosed when one of said evaporator conductance-area product and 
said evaporator heat exchanger approach exceeds said predetermined threshold value 
for said CQ. 

1 1 . The fault detection system according to claim 6 wherein condenser tube 
fouling is diagnosed when one of said condenser conductance-area product and said 
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condenser heat exchanger approach exceeds said predetermined threshold value for 
said CQ. 

12. The fault detection system according to claim 6 wherein an internal 
compressor internal fault is diagnosed when one of said compressor coefficient of 
performance, said isentropic efficiency, and compressor power draw exceeds said 
predetermined threshold value for said CQ. 

13. The fault detection system according to claim 6 wherein a motor 
transmission fault is diagnosed when one of said overall coefficient of performance, 
said motor/transmission efficiency, and motor power draw exceeds said 
predetermined threshold value for said CQ. 

14. The fault detection system according to claim 6 wherein said 
predetermined threshold includes upper and lower critical value levels for each said 
characteristic quantity, said upper and lower critical value levels being determined in 
accordance with a rated precision error of a sensor used to acquire a selected said 
measured input used to calculate said characteristic quantities. 

15. A method for detecting faults in a chiller subsystem of a facility cooling 
system, comprising: 
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3 providing plural base data sets of measured inputs, each said set of measured 

4 inputs including sensor data for a different base operating condition of the chiller; 

5 computing a set of characteristic quantities for each of said plural base data 

6 sets; 

7 storing said plural sets of characteristic quantities in a memory means; 

8 providing a test data set for a test operating condition of the chiller; 

9 computing a set of characteristic quantities from said test data set; 

! o producing an interpolated base set of characteristic quantities from among said 

1 1 plural sets of characteristic quantities stored in said memory means using said test 

12 data set; 

13 detecting a fault when at least one characteristic quantity exceeds a 

1 4 predetermined threshold range; 

! 5 diagnosing a fault in relation to which ones of said characteristic quantities 

1 g exceed said predetermined threshold range. 

1 16. The method according to claim 15 wherein said predetermined 

2 threshold range includes upper and lower critical value levels for each said 

3 characteristic quantity, said upper and lower critical value levels being determined in 

4 accordance with a rated precision error of a sensor used to acquire data items 

5 contained in said test data set. 
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17. The method according to claim 15 wherein said measured inputs are 
selected from the group comprising: {chilled water supply temperature; chilled water 
return temperature; condenser water supply temperature; chilled water flow rate; 
condenser water flow rate; evaporator saturation temperature; condenser saturation 
temperature; and compressor discharge temperature. 

18. The method according to claim 1 5 wherein said characteristic quantities 
are selected from the group comprising: evaporator heat exchanger approach; 
condenser heat exchanger approach; chilled water temperature difference; condenser 
water temperature; condenser conductance-area product; condenser conductance-area 
product; isentropic efficiency; motor/transmission efficiency; overall coefficient of 
performance; and compressor coefficient of performance. 

19. The method according to claim 1 8 wherein condenser water flow rate 
abnormalities are diagnosed when one of said condenser conductance-area product, 
said condenser heat exchanger approach, and said condenser water temperature 
difference exceeds said predetermined threshold value for said CQ. 

20. The method according to claim 18 wherein chilled water flow rate 
abnormalities are diagnosed when one of said evaporator conductance-area product, 
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3 said evaporator heat exchanger approach, and said chilled water temperature 

4 difference exceeds said predetermined threshold value for said CQ. 

1 21. The method according to claim 18 wherein water flow rate 

2 abnormalities in said chilled water flow are detected when one of said evaporator 

3 conductance-area product, said evaporator heat exchanger approach, and said chilled 

4 water temperature difference exceeds said predetermined threshold value for said CQ. 

1 22. The method according to claim 1 8 wherein evaporator tube fouling is 

2 diagnosed when one of said evaporator conductance-area product and said evaporator 

3 heat exchanger approach exceeds said predetermined threshold value for said CQ. 

1 23. The method according to claim 18 wherein condenser tube fouling is 

2 diagnosed when one of said condenser conductance-area product and said condenser 

3 heat exchanger approach exceeds said predetermined threshold value for said CQ. 

1 24. The method according to claim 18 wherein an internal compressor 

2 internal fault is diagnosed when one of said compressor coefficient of performance, 

3 said isentropic efficiency, and compressor power draw exceeds said predetermined 

4 threshold value for said CQ. 
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25. The method according to claim 1 8 wherein a motor transmission fault 
is diagnosed when one of said overall coefficient of performance, said 
motor/transmission efficiency, and motor power draw exceeds said predetermined 
threshold value for said CQ. 

26. A fault detection and diagnosis system for a HVAC system including 
at least two chillers, each chiller being equipped with sensors for measuring the 
performance of one or more condenser, a compressor, an evaporator, and a chiller, 
said fault detection system comprising: 

processing means for calculating characteristic quantities (CQ's) from a 
plurality of measured inputs for each of the chillers, one chiller being designated as 
a base-case chiller; 

means for detecting a fault when a difference between CQ values for the base 
case chiller and CQ values for other ones of the chillers exceeds a predetermined 
threshold range for at least one of said CQ values; 

a fault classifier for classifying a detected fault based on which ones of said 
actual CQ values exceed said threshold range. 

27. The fault detection system according to claim 26 wherein said CQ's are 
selected from the group comprising: evaporator heat exchanger approach; condenser 

"heat exchanger approach; chilled water temperature difference; condenser water 
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temperature; condenser conductance-area product; condenser conductance-area 
product; isentropic efficiency; motor/transmission efficiency; overall coefficient of 
performance; and compressor coefficient of performance. 

28. The fault detection system according to claim 27 wherein condenser 
water flow rate abnormalities are diagnosed when one of said condenser conductance- 
area product, said condenser heat exchanger approach, and said condenser water 
temperature difference exceeds said predetermined threshold value for said CQ. 

29. The fault detection system according to claim 27 wherein chilled water 
flow rate abnormalities are diagnosed when one of said evaporator conductance-area 
product, said evaporator heat exchanger approach, and said chilled water temperature 
difference exceeds said predetermined threshold value for said CQ. 

30. The fault detection system according to claim 27 wherein chilled water 
flow rate abnormalities are diagnosed when one of said evaporator conductance-area 
product, said evaporator heat exchanger approach, and said chilled water temperature 
difference exceeds said predetermined threshold value for said CQ. 

3 1 . The fault detection system according to claim 27 wherein evaporator 
tube fouling is diagnosed when one of said evaporator conductance-area product and 
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3 said evaporator heat exchanger approach exceeds said predetermined threshold value 

4 for said CQ. 

1 32. The fault detection system according to claim 27 wherein condenser 

2 tube fouling is diagnosed when one of said condenser conductance-area product and 

3 . said condenser heat exchanger approach exceeds said predetermined threshold value 

4 for said CQ. 

1 33. The fault detection system according to claim 27 wherein an internal 

2 compressor internal fault is diagnosed when one of said compressor coefficient of 

3 performance, said isentropic efficiency, and compressor power draw exceeds said 

4 predetermined threshold value for said CQ. 

1 34. The fault detection system according to claim 27 wherein a motor 

2 transmission fault is diagnosed when one of said overall coefficient of performance, 

3 said motor/transmission efficiency, and motor power draw exceeds said 

4 predetermined threshold value for said CQ. 

1 35. The method according to claim 27 wherein said predetermined 

2 threshold includes upper and lower critical value levels for each said characteristic 

3 quantity, said upper and lower critical value levels being determined in accordance 
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4 with a rated precision error of a sensor used to acquire a selected said measured input 

5 used to calculate said characteristic quantities. 

1 36. The fault detection system according to claim 26 further comprising: 

2 means for producing interpolated base-case CQ values from base-case chiller 

3 for a given set of measured inputs; 

4 wherein said fault detecting means detects a fault when a difference between 

5 said interpolated base-case CQ values and CQ values for other ones of the chillers 

6 exceeds a predetermined threshold range for at least one of said CQ values. 
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FIG. 1 
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FIG. 3 

Procedure LMTDe_(T_chws,T_chwr,T_evap:LMTDe,T_we,ln_arg_e) 
T_we = (T_chws +T_chwr)/2 

Inarge = max(0,(T_chwr-T_evap)/(T_chws-T_evap)) 
If ((T_chws = T_chwr) OR (ln_arg_e=0)) Then 
LMTDe = T_we - T_evap 
Else 

If (T_evap > T_chws) Then Call Error 
If (T_chws > T_chwr) Then Call Error 

LMTDe_real = ((T_chwr-T_evap)-(T_chws-T_evap))/ln(ln_arg_e) 
LMTDe = max(0,LMTDc_real) 

Endlf 
End 



FIG. 4 

Procedure LMTD_c(T_cws,T_cwr,T_cond:LMTDc,T_wc,ln_arg_c) 
T_wc = (T_cws + T_cwr)/2 

In_arg_c - max(0,(T_cond-T_cws)/(T_cond-T_cwr)) 
If (T_cwr = T_cws) OR (ln_arg_c =0) Then 
LMTDc = T_cond - T_wc 
Else 

If (Tcond < T cwr) Then Call Error 
If (T_cwr < T_cws) Then Call Error 

LMTDc_real = ((T_cond-T_cws)-(T_cond-T_cwr)) / ln(ln_arg_c) 
LMTDc = max(0,LMTDc_real) 

Endlf 
End 
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FIG. 5 

time = tablerun# 

T_cws = lookup(time, #T_cws) "F" 
Gpmcw - lookup(time, #Gpmcw) "gpm" 
T_cond = lookup(time, #T_cond) "F" 

FIG. 6 

T_chws = lookup(time, #T_chws) "F" 
T_chwr = lookup(time, #T_chwr)"F" 
Gprachw = lookup(time, #Gpmchw) "gpm" 
T_evap = lookup(time, #T_evap) "F" 

FIG. 7 

T_2 = lookup(time, #T_2) "F" 
Power = lookup(time, #Power) "kW" 



FIG. 8 

Q_evap = Gpmchw * ft3hr\gpm * dens_chw * Cpw * (T_chwr - T_chws)/1 2000 "ton" 
Q_evap = UAe * LMTDe /12000 "ton" 

Q_evap = Gpm R * ft3hr\gpm * dens_R* (h_l - h_4)/12000 "ton" 

h_l = Enthalpy(R22, x = 1 , P = P_evap) "Btu/Ibm" 

v_l = Volume(R22, x= 1, T- T_evap) "ft A 3/lbm" 

T_evap = Temperature(R22,x = 0,P = P_evap) "F" 

APPRe = T_chws - T evap "F" 

CHWTD = T_chwr - T chws "F" 
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FIG. 9 

P_comp = Gpm_R * ft3hr\gpm * dens_R * (h_2 - h_l)/120Q0 * kW\ion "kW" 
P_comp = P_ideal/eta "kW 

P_ideal - Gpm_R * ft3hr\gpm * dens_R * (h_2s - h_l)/12000 *kW\ton "kW" 
s_l = Entropy(R22, x - 1, P = P_evap) "Btu/lbm-R" 
s_2s = s_l "Btu/lbm-R" 

h_2s = Enthalpy(R22, s = s_2s, P = P_cond) "Btu/lbm" 
COP_overall = Q_evap /(Power / kW\ton) 
COP_comp = Q_evap /(P_comp / kW\ton) 
COP_ideal = Q_evap /(P_ideal / kW\ton) 
EFFMOT = P_comp/Powcr 
P_ratio = P_cond/P_evap 



C& 02344908 2001-04-23 



FIG. 10 

"CONDENSER EQUATIONS:" 

Q_cond = Gpmcw * ft3hr\gpm * dens cw * Cpw * (T_cwr - T_cws)/12000 "ton" 
Qcond = UAc * LMTDc /12000 "ton" 

Q_cond = Gpm_R * ft3hr\gpm * dens_R * (h_2 - h_3)/12000 "ton" 

h_3 = Enthalpy(R22, x = 0, P = P_cond) "Btu/lbm" 

v_3 = Volume(R22, x = 0, P = P_cond) "ft A 3/lbm" 

T_cond = Temperature(R22, x = 0, P = P_cond) "F" 

T_2 = Temperature(R22, h = h_2, P = P_cond) "F" 

v_2 = Volume(R22, h - h_2, P = P_cond) "ft A 3/lbm" 

APPRc = T cond - T_cwr "F" 

CWTD = T_cwr - T_cws "F" 

FIG. 11 

"EXPANSION VALVE EQUATION:" 
h_4 = h_3 "Btu/lbm" 



FIG. 12 

"PROPERTY & CONVERSION DATA:" 

dens_chw = Density(Water, T=(T_chws+T_chwr)/2, P=14.7) "lbm/ft A 3" 

dens_cw = Density(Water, T=(T_cws+T_cwr)/2, P=14.7) "lbm/ft A 3" 

dens_R = Density(R22, T=(T_cond+T_evap)/2, P=(P_cond +P_evap)/2) "lbm/ft' 

Cpw= 1.0 "Btu/lbm-R" 

ft3hr\gpm = eonvert(gal,ft3)/convert(min,hr) 

kW\ton = convert(ton,kW) 
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FIG. 13 
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FIG. 15 




Input layer Hidden layer Output layer 
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FIG. 18 




FIG. 19 
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FIG. 21 
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FIG. 22 




-15 -10 -5 0 5 

Isentropic Efficiency Deviation [%] 



CA 02344908 2001-04-23 



FIG. 23 
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FIG. 26 
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chiller. The system has a thermodynamic preprocessor for 
calculating characteristic quantities (CQ's) from a plurality of 
measured inputs. Also provided is a base-case lookup table for 
storing plural first sets of CQ values generated by the 
thermodynamic preprocessor in an initial period in which fault-free 
operation of the HVAC system is assumed, a first set of CQ values 
being generated for each of plural different measured input values. 



DESCRIPTION - An interpolator interpolates a set of base-case CQ 
values from the first sets of CQ values stored in the base-case 
lookup table for a given set of measured inputs. A fault detector 
detects a fault when a difference between actual CQ values and base- 
case CQ values exceeds a predetermined threshold value for at least 
one of the CQ values. A fault classifier classifies a detected 
fault based on which ones of the actual CQ values exceed the 
interpolated CQ values. The processing unit uses measured inputs 
selected from the group comprising: chilled water supply 
temperature; chilled water return temperature; condenser water 
supply temperature; chilled water flow rate; condenser water flow 
rate; evaporator saturation temperature; condenser saturation 
temperature; and compressor discharge temperature. The interpolator 
may be a General Regression Neural Network, a Probabilistic Neural 
Network, or a Back Propagation Network. 
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